Charge fluctuation induced dephasing in a gated mesoscopic interferometer 
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The reduction of the amphtude of Aharonov-Bohm osciUations in a ballistic one-channel meso- 
scopic interferometer due to charge fluctuations is investigated. In the arrangement considered the 
interferometer has four terminals and is coupled to macroscopic metallic side-gates. The Aharonov- 
Bohm oscillation amplitude is calculated as a function of temperature and the strength of coupling 
between the ring and the side-gates. The resulting dephasing rate is linear in temperature in 
agreement with recent experiments. Our derivation emphasizes the relationship between dephasing, 
ac-transport and charge fluctuations. 



I. INTRODUCTION 

Dephasing processes suppress quantum mechanical in- 
terference effects and generate the transition from a mi- 
croscopic quantum coherent world in which interference 
is crucial to a macroscopic world characterized by the 
absence of (quantum) interference effects. Mesoscopic 
systems are neither entirely microscopic nor macroscopic 
but at the borderline between the two. Clearly, there- 
fore, dephasiip|E-,processes play a central role in meso- 
scopic physicsErla. At low temperatures, it is thought that 
the predominant process which.|generates dephasing, are 
electron-electron interactionsd'aB. In this work, we inves- 
tigate a ballistic Aharonov-Bohm (AB) interferometer, in 
which electrons (in the absence of interactions) are sub- 
ject only to forward scattering processes (see Fig. (|l|)). 

Our work is motivated by the following questions: 
A mesoscopic conductor connects two or more electron 
reservoirs: inside a reservoir screening is effective and 
electron interactions are of little importance. In con- 
trast, inside the mesoscopic structure screening is poor 
and interactions are important. Thus the process of a 
carrier entering or leaving the mesoscopic structure is 
essential. We ask how this process affects the dephas- 
ing. In standard treatments of dephasing the conductor 
is considered to be charge neutral and the elementary 
excitations are electron-hole scattering processes. In a 
finite size mesoscopic conductor, we can, however, have 
a hole in a reservoir and only an additional electron in 
the conductor, or we can have an electron inside the con- 
ductor and a hole on a nearby capacitor (see Fig. |l]) . As 
a consequence the conductor is charge neutral only when 
its surroundings are taken into account (reservoir banks, 
nearby capacitors). 

A second question we seek to answer is the following: 
instead of calculating a dephasing rate it is desirable to 
find a way to directly evaluate the quantity of interest 
(here the conductance) . In small mesoscopic systeros the 
dephasing rate might be a sample specific quantityO and 
there would be little justification in using an ensemble av- 
eraged dephasing rate even if we are interested only in the 
ensemble averaged conductance. Clearly to answer such 
conceptual questions it is useful to have a model which is 



as simple as possible. In this paper we theoretically in- 
vestigate dephasing of AB-oscillations in a ballistic ring 
with a single transport channel. The one-channel limit 
is of actual experimental interest (see e.g. Refs. |7|-||). 
Our idealized setup consists of an AB-ring with four ter- 
minals, the arms of the ring being capacitively coupled 
to lateral gates (see Fig. (|])). For a recent experiment 
on a ballistic (two-terminal) AB-ring with lateral gates 
coupled to both arms see Ref. |lO|. 
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FIG. 1. The figure shows the four-terminal AB-ring 
threaded by a magnetic fiux. The two arms of the ring are 
each coupled to a side-gate via a capacitance Co (G = A,B). 
We consider junctions which are perfectly transmitting and 
divide the incoming current into the upper and lower branches 
of the ring. The system then is the electronic equivalent of an 
optical Mach-Zehnder interferometer (MZI) . The total charge 
in a Gauss sphere 57 drawn around the system of gates and 
ring is assumed to be zero, implying that current in the system 
is conserved. It is assumed that each arm is characterized by 
a single potential Ua (or Ub) and the charge -^Qa (or +Qb)- 

The structure we examine has no closed orbits: As a 
consequence in its equilibrium state it exhibits no persis- 
tent current and in the transport state there is no weak 
localization correction to the conduietance. It exhibits, 
however, an Aharonov-Bohm effecttll due to superpo- 
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sition of partial waves in the out-going final quantum 
channel. In fact this is the situatiop, discussed in the 
original work of Aharonov and BohrnO. It is also some- 
times assumed in mesoscopic physics without a detailed 
specification of the conditions (multi-terminal geometry, 
absence of backscattering) which are necessary for inter- 
ference to appear only in the final outgoing channel. The 
system investigated here is the electric analog of an opti- 
cal interferometer in which the path is divided by forward 
scattering only. An example of such an arrangement is 
the Mach-Zehnder interferometer (MZI). In the MZI the 
sample specific AB oscillations are a consequence of su- 
perposition in the outgoing scattering channel only. We 
calculate the effect of internal potential fluctuations on 
the linear response dc-conductance in the ring as a func- 
tion of the strength of the coupling between ring and 
gates and of temperature. In this approach the conduc- 
tance and the dephasing time are not calculated sepa- 
rately but a dephasing time will appear in the expression 
for the conductance in a natural way. It quantifies the 
degree of the attenuation of AB-oscillations due to ran- 
domization of the phases of the electrons going through 
the ring (as opposed to attenuation due to thermal av- 
eraging). From our calculations we find that the coher- 
ent part of the conductance is diminished by a factor 
1 — t/t^ « exp(— t/t^) relative to the ideal case due to 
temperature and coupling to the gates. Here r is the 
traversal time for going through one arm of the ring. In 
the temperature regime ^ kT ^ AE {AE is the sub- 
band threshold) addressed in our calculations we find a 
dephasing rate t^^ linear in temperature. In an experi- 
ment on a two-terminal AB-ring, a dephasing rate linear 
temperature was recently measured by Hansen et al. 



In our model the dephasing is due to inelastic scatter- 
ing of electrons from charge fluctuations in the arms of 
the ring. We treat the gates as macroscopic entities with 
perfect screening. The carrier dynamics in the gates is 
irrelevant for the discussion presented here. They do not 
represent an external bath or dephasing agent. The ir- 
reversible source necessary for dephasing is given by the 
electron dynamics of the ring itself: the phase and en- 
ergy of a carrier exiting into a contact are unrelated to 
the phase and energy of a carrier entering the sample. 

Our approach is similar in spirit to the one used in 
Ref. |l^ where dephasing due to charge fluctuations was 
discussed for two coupled mesoscopic structures. To sim- 
plify the discussion we assume that it is possible to draw 
a Gauss sphere around the system of gates and ring such 
that all electric field lines emanating within the sphere 
also end in it (see Fig. (|l])). This implies that the total 
charge in the sphere is zero at any time. However, it is 
possible to charge up one part of the system (an arm of 
the ring) relative to another part (the nearby gate) cre- 
ating charge dipoles. Charge fluctuations in the arms of 
the ring lead to fluctuations of the effective internal po- 
tentials. Electrons going through the ring are exposed to 
these potential fluctuations and scatter inelastically. The 



strength of the coupling between the gates and the arms 
determines the amount of screening and thus the strength 
of the effective electron-electron interaction. When the 
capacitance C between arm and gate becomes very large, 
the Coulomb energy e^/C of the system goes to zero and 
gate and arm are said to be decoupled. 

The main goal of this work is the calculation of the 
dc-conductance of the ring taking into account the ef- 
fects of the gate-mediated interactions. Applying the ac- 
scattering approach of Ref. |l3| we start by calculating 
the dynamic conductance matrix of our system. From 
the real (dissipative) part of the dynamic conductance 
matrix element relating the current in one of the gates 
to the voltage applied to the same gate we can find the 
spectrum of the equilibrium potential fluctuations in the 
nearby arm of the ring via the fluctuation dissipation 
theorem. When calculating the dc-conductance we sta- 
tistically average over the scattering potential assuming 
that the potential has a vanishing statistical mean and 
that the spectrum of fluctuations is given by the spec- 
trum of equilibrium fluctuations. Taking into account 
interactions in the manner outlined above results in an 
attenuation of the amplitude of AB-oscillations of the 
statistically averaged conductance. 

In the next section we derive the scattering matrix for 
the MZI. Inelastic scattering from internal potential fluc- 
tuations is taken into account. In Section IV we deter- 



mine the internal potential distribution of the ring and 
then go on to calculate the admittance matrix. We will 
show that by taking into account screening effects a cur- 
rent conserving theory for the system consisting of the 
ring and the two gates can be formulated. In the fol- 
lowing (Section ^ we calculate the dc-conductance and 
investigate the influence of equilibrium fluctuations on 
dc-transport. 



II. THE MACH-ZEHNDER INTERFEROMETER 

We consider an MZI with a single transport channel. 
An electron arriving at one of the two intersections com- 
ing from a reservoir can enter either of the two arms 
of the ring, but can not be reflected back to a reser- 
voir. An electron coming to the intersection from the 
ring will enter one of the reservoirs. The amplitudes for 
going straight through the intersection and for being de- 
flected to the adjacent lead in the forward direction are 
t ~ \fT and r = i-y/i? respectively, where T + R — 1. 
Transmission through the intersections is taken to be in- 
dependent of energy. In the remainder, we assume sym- 
metric intersections, that is i? = T= 1/2 in (Eq. |^). 
Due to the potential fluctuations in the arms of the ring 
a carrier can gain or loose energy. This process is de- 
scribed by a scattering matrix Sg{E',E) for each arm 
which depends on both the energy E of the incoming 
and the energy E' of the exiting carrier. The scatter- 
ing matrix Sg{E\E) thus connects current amplitudes 
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at a junction of the ring incident on the branch to the 
ampHtudes of current at the other junction leaving the 
branch. We have a matrix Sg{E',E) for the upper arm 
{Sa{E\E)) and one for the lower arm {Sb{E',E)). As 
a consequence of the inelastic transitions in the arms of 
the ring the full scattering matrix Sa/iiE' , E) describ- 
ing transmission through the entire interferometer from 
contact f3 to contact a is also a function of two energy 
arguments. This scattering matrix can be found by com- 
bining the scattering matrices for the two arms with the 
amplitudes for going through the intersections following 
a specified path. Due to the geometry of the system we 
have the symmetries Si3^+b{E' , E) = S24,+b{E' , E) and 
S32,+b{E' , E) = Sii^+siE' , E). In addition the scatter- 
ing matrix elements calculated for the system in a mag- 
netic field B are related to the matrix elements found at 
an inverted field -B by Sai3,+B{E', E) = Spa,-B{E', E). 
All elements of the scattering matrix can then be found 
from the three elements given below: 

Si:i{E',E)^S2i{E\E)^iVTR (1) 

X {SA{E\E)e-'-^^ + Sb{E' ,E)e+''^^) , 
Su{E',E) = -RSA{E',E)e-"^'^ +TSB{E',E)e+"^", 
S2s{E',E) =TSA{E',E)e-"^^ - RSb{E' , E)e+"^'' . 

Here $g is the magnetic phase picked up by a particle go- 
ing through arm G clockwise. Then ^a + ^b ~ 27r$/$Oi 
where $ the flux through the ring and <f>o is the fiux quan- 
tum. All scattering matrix elements Sa/3{E' , E) with 
|a — /3| < 1 are zero since transport through the junc- 
tions takes place only in forward direction. The on-shell 
(one-energy) scattering matrix elements for the (free) 
one-channel interferometer in the absence of gates which 
we denote by Sf^ {E) are found by replacing Sg{E' , E) in 

Eq. 1) by Sg{E',E) = 5{E'~E)s'§\e) with S^^\e) = 
explikELc)- Our first task is now to determine the scat- 
tering matrices Sg^E' , E) for the arms of the ring. 



III. S-MATRIX FOR A TIME-DEPENDENT 
POTENTIAL 

In this section we calculate the scattering matrix 
for the interacting ring system. We first solve the 
Schrodinger equation for a single branch of the interfer- 
ometer using a WKB approach. The amplitude for a 
transition from energy E to energy E + Tilj of a particle 
passing through this arm and the corresponding scatter- 
ing matrix element are determined. Subsequently the 
scattering matrices for the two arms are included into a 
scattering matrix for the full interferometer. A WKB ap- 
proach similar to ours has been used previously to discr^ 
photon-assisted transport in a quantum point contactEfl 
or to the investigation of traversal times for tunnelingll3. 
The influence of a time-dependent bosonic environment 
on transport through a QPC was addressed in Refs. |l^ 
and |l^ also applying a WKB ansatz. 



The gate situated opposite to arm G {G = A, B) is 
assumed to be extended over the whole length of this 
arm. Fluctuations of the charge in the gate capacitively 
couple to the charge in the neighboring arm of the MZI 
and influence electron transport through this arm. This 
interaction effect is taken care of by introducing a time- 
dependent potential VG{x,t) into the Hamiltonian 



H = - 



2m* dx^ 



+ EG + VGix,t), 



(2) 



for arm G. Here Eg is the sub-band energy due to the 
lateral confining potential of the arm and m* is the ef- 
fective mass of the electron. We make the assumption 
that the fluctuating potential factorizes in a space- and 
a time-dependent part, writing Vg{x, t) = hG{x) e UG^t)- 
For the ballistic structure considered here the internal 
potential is a slowly varying function of x. For practical 
calculations we will however often employ a rectangular 
potential barrier, (JiGix) = const if < a; < Lg, where 
Lg is the length of arm G) . Using a spacejndependent 
internal potential is a valid approximationllj at least in 
the low frequency limit ujtg ^ 1 where a passing electron 
sees a constant or slowly changing barrier. We have here 
introduced the traversal time tg = Lg/vg,f where vg,f 
is the Fermi velocity in arm G. We will show in Section 
how the potentials 



UG{t) 



duj 
2^ 



UG{uj)e 



-iujt 



(3) 



and their spectra can be determined in a self-consistent 
way. To solve the Schrodinger equation with the Hamil- 
tonian Eq. (^ we make the ansatz 



(4) 



where fc^ ^ = y/2m* [E — Eg)/Tt, and rG{x,t) is the ac- 
tion due to inelastic scattering. We will omit the index G 
of the wave vector (or of the velocity ve,g = ftkG.E/fn*) 
from now on and only write it when the distinction be- 
tween wave vectors in different arms is important. It is 
assumed that transmission is perfect: the potential fluc- 
tuations cause only forward scattering. This is justified 
if all energies relating to the fluctuating potential are 
much smaller than the Fermi energy Ep, in particular 
Hlu < Ep- 

In determining rG{x,t) we will take into account cor- 
rections up to the second order in the potential. We write 

rG{x, t) = rG,i{x, t) + TG^ix, i), where 



rGAx.t) 



duj 
— ( 
27r 



'rG,i(a;, w) 



(5) 



is linear in the perturbing potential and 
duji f doj 



rG,2{x,t) 



2n 



2-n 



ie-»(-i+"2)*rG,2(a;,c^i,c^2) (6) 
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is a second order correction. The linear term was cal- 
culated in Ref. |l5| for a general form of hc{x). The 
corresponding expression for the term quadratic in the 
potential is readily found but is quite cumbersome. We 
here give ro,i{x,uj) and rG,2{x,LOi,uj2) for the case were 
hcix) is a rectangular barrier of length Lq: 



(7) 



2m*Vp 

Here re,! gives the contribution to the action due to ab- 
sorption or emission of single modulation quantum huj, 
while r(3_2 corresponds to the absorption and emission of 
two modulation quanta huji and huj2- We now proceed 
to the formulation of the scattering problem in terms of 
a scattering matrix with elements of the form Sg{E' , E) 
describing transitions between states at different ener- 
gies. The amplitude tQ{E',E) for a transition from a 
state with energy to a state with energy E' of an elec- 
tron is found from the boundary condition sX x = Lq, 
'^%{LG,t) = Xe {^G,t) -For the matching we expand the 
WKB wavefunction (see Eq. (^)) in a; = Lq to the second 
order in the perturbing potential 



^-ilx,*) =e 



l + j^rcALct) (8) 
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Furthermore the wavefunction at the right-hand side of 
the barrier (outside the fluctuating potential region) is 



iE't/h 



(9) 



In principle, also the derivatives of the wavefunctions 
should be matched. Here we describe transmission 
through the fluctuating potential region as reflectionless 
which is accurate up to corrections of the order of hco/Ep. 
To determine the transmitted wave with the same accu- 
racy it is sufficient to match amplitudes only. The trans- 
mission amplitude is found by Fourier transforming the 
WKB wavefunction Eq. (H) and comparing with Eq. (|^). 
The transmission amplitude can be expressed in terms of 
the phase raiLc, t). To second order in the potential we 
have 



tG{E',E) = e 



2TThSie) + -rGALG,e/h) 



duj 
2^ 



TrG,2iLG,UJ,s/h~Uj) (10) 

n 



1 

2^^ 



rG,iiLG,uj)rGs iLG,s/h-uj) 



where e = E' — E. The scattering matrix connecting 
incoming wave amplitudes (at a; = 0) to outgoing wave 



amplitudes (at x = Lg) is related to the transmission 
amplitude tG{E',E) through 



SgIE', E) = e'^'^'^''tG[E\ E) 



(11) 



While the transmission amplitude tG{E',E) was deter- 
mined through the continuity of the wavefunction in the 
point X = Lg (tG{E' , E) thus connects amplitudes at the 
same point), the scattering matrix Sg{E',E) connects 
amplitudes in a; = to amplitudes in x = Lg- This dif- 
ference in the definitions of the two quantities leads to 
the phase factor ex.p{ikE' Lg) in Eq. (^. The scattering 
matrix as it is derived here a priori relates wavefunction 
amplitudes and not current amplitudes. To be consistent 
with the usual definition of the scattering matrix as a re- 
lation between current amplitudes, Sg{E',E) should be 
multiplied by \/ ve' /ve- This factor, however, is of the 
order Tuju/Ef and can thus be neglected. The scattering 
matrices found to describe a single arm can now be in- 
tegrated into the full scattering matrix for the MZI (see 
Eq. (0)). 



a (E+liCi)) 
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h(x)eU(t) 



FIG. 2. Scattering states *g and ^b^r^ (see Eq. (^) 
with amplitudes a{E) and a{E±Tiui) respectively due to elec- 
trons incident from the left at energies E and E ± Tiuj are 
indicated in the figure. For didactic purposes the special case 
of a harmonically oscillating barrier U{t) cx cositut) is con- 
sidered and only first order side-bands are drawn. In the 
rest of the paper we discuss the case of a randomly oscillat- 
ing barrier and include second order corrections. The scat- 
tering state which may be described by a simple (outgoing) 
plane wave at energy E with amplitude h{E) to the right 
of the barrier is emphasized in the figure. The amplitudes 
a{E) and b{E) are related through the scattering matrix via 
KE) = Ea=o,± ^(E, E + ahuj)a{E + ahuj). 

In the discussion presented here the transmission of 
the carrier through the fluctuating potential region is de- 
scribed as a unitary scattering process. The "final" scat- 
tering channels are always open. We emphasize that up 
to now we have investigated a perfectly coherent process. 
Decoherence in our model will be introduced through the 
statistical averaging (cf. Sec. 0). Our next task is to 
find the statistical properties of the potential fiuctua- 
tions. These fiuctuations can be found from the dynamic 
conductance matrix via the fluctuation-dissipation theo- 
rem. 
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IV. POTENTIAL FLUCTUATIONS 



In this section we proceed to the calculation of the 
admittance matrix Gik{oj) = dli^uj / dVu^u for the joint 
system of interferometer and gates. We concentrate on 
the hmit hoj <ti kT <ti Ep. The dynamic conductance 
matrix Gik{oj) is a 6 x 6-matrix (/, fc = 1,2,3,4,^,5), 
Ii^i^ and Vk^u denoting respectively the current measured 
at and the voltage applied to one of the four contacts of 
the ring or to one of the two gates. We use the follow- 
ing convention for the indices: Lower case Roman indices 
can take the values 1, 2, 3, 4, A, B, the Greek indices a, (3 
take the values 1,2,3,4 while upper case Roman indices 
G, H are A, B. We will first calculate the matrix elements 
Gcci'^) from which, via the fluctuation dissipation the- 
orem we can derive the spectra of potential fluctuations 
in the two arms. These will later be needed in the discus- 
sion of the decoherence of AB-oscillations. The remain- 
ing elements of the conductance matrix and the resulting 
total conductance matrix are given in Appendix ^ The 
elements of the conductance matrix obey the sum-rules 
J2iGiki(^) = and J2k'^ik{^) = reflecting gauge in- 
variance and current conservation respectively. A prob- 
lem closely related to the one addressed here is concerned 
with the calculation of ac-transport properties of a bal- 
listic wire-attached to reservoirs and capacitively coupled 



to a gateEfl. Contrary to the classical calculations done 
for a wire in Ref. ^ the ac-scattering approach allows us 
to take into account the quantum nature of the system 
investigated here as manifested in the AB-oscillations. 

The ac-properties and the potential distribution which 
are of interest here depend not only on the mesoscopic 
conductor but also on the properties of the external cir- 
cuit. Here we consider the case where all external current 
loops exhibit zero impedance. This requires some expla- 
nation since especially voltages at the gates are typically 
controlled with the help of an external impedance. How- 
ever, what counts in our problem is the range of fre- 
quencies up to the traversal time, whereas the external 
impedance might be very large only in a very narrow fre- 
quency range around lo — 0. Thus we are justifled to 
consider in the following a zero-impedance external cir- 
cuit. 

In order to obtain the conductance matrix from an 
ac scattering approach we need the effective internal po- 
tential eUcit) in arm G. The internal potential eUcit) 
is related to the total charge Qoit) in the same region 
through Qoit) = CciUait) - Vcit)) where Vb(i) is the 
voltage applied to gate G and Cg is a geometrical capaci- 
tance characterizing the strength of the coupling between 
arm and gate. The total charge Qait) consists of a con- 
tribution due to injection from the contacts labeled QqII) 
and a screening part Q^it), thus Qait) = QQ{t) — QQ{t). 
We will now assume that a voltage Va{t) is applied to 
contact a while Vf3{t) = for a 7^ /3. 

First we consider the charge densityli3 injected into the 
arm G due to a modulation of the voltage at contact a 



assuming a fixed internal potential Ug- The charge dis- 
tribution in the sample can be expressed through the 
Fermi-fleld 

Hr,t)^Yl I -^^e-^'''/'^Mr;E)aUE), (12) 

which annihilates an electron at point r and time t. Here 
(r; E') is a scattering state describing carriers with en- 
ergy E incident from contact a. The charge density in 
the ring at point r and time t is p(r, t) — e^^(r, i)^'(r, t). 
Fourier transforming with regard to time and quantum 
averaging we get paij^uj) = {pa{r,Lu))^ where 



P„(r,.) = e5:|-— 



dE 



(13) 



X E)^l^p{Y- E + hLo){al{E)af3{E + hcu)). 

The average charge may be split into an equilibrium part 
p'^"^ (r, Lu) and a contribution due to the external voltage 
Spa{r,uj): 



p{r,uj) = p^°\r,uj) + Spair,uj). 



(14) 



When calculating the quantum average of the charge den- 
sity operator the effect of the external voltage is taken 
into account through the modified distribution function 
for charge carriers coming in from reservoir a. The distri- 
bution for contact a to linear order in the applied voltage 



{dl{E)a„{E + hu)) 



(15) 



:d{hLu)U{E) + -V^,^F{E,Lu), 
n 



where Va^^ is the Fourier component to frequency lo of 
the voltage Va{t) and F{E,u!) is defined through 



F{E,Lu) = 



U{E)- U{E + hLo) 

huj 



(16) 



Carriers in the other reservoirs are Fermi distributed 
{{ai{E)af3iE + hu)) = S{hu;)S^i3f^{E) for a ^ ^ or 
/? ^ fi). The scattering states ipa{ic;E) in the arms of 
the interferometer for a constant internal potential are of 
the form ipa{r;E) = ^Qx(r^) exp(zA:£;a;-|-i<I>G(a::)), where 
Aa = i\/R or Aa = VT depending on the arm and the 
injecting contact (cf. Eq. (|l])). As in most of the paper 
we will in the following use R = T = 1/2. Furthermore 
$g(x) is the magnetic phase acquired going through arm 
G to point X and x(rx) is the transverse part of the 
wavefunction. The simple form of the scattering states 
in the arms is a consequence of the absence of backscat- 
tering in the intersections. The injected charge 5pa{x, uj) 
is the part of the total charge Eq. ( |l^ ) proportional to the 
non-eq uilibrium contribution to the distribution function 
Eq. (Ila). Substituting the expressions for the scattering 
states into Eq. (p^), using Eq. (|l^) and integrating over 
r_L we find 
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e 



dE 



(17) 



where we have used |^aP = 1/2. To find the total charge 
Qg a(^) injected into arm G of the MZI wc integrate over 

the length of the arm Qq ^{ut) 
forming the integration we get 



dxSpa{x,Lu). Per- 



lG,a\ 



H = ^ / dEFiE^u) ( - ) (1 



(18) 



In the limit fko/kT < 1 we have fdEF{E,uj) « 1. We 
can rewrite the charge as Qg a('^) ~ '^'^i^Ga{^)Va,uj where 
we have introduced the injectivity i'gq(^)j defined as 



(19) 



Here tg = Lg/vp is the traversal time through arm G. 
Now if interactions are taken into account, the excess in- 
jected charge will induce a shift in the effective internal 
potential, which in turn gives rise to a screening charge. 
This screening charge is proportional to the internal po- 
tential euG{^) and to the total charge density available 
for screening vg{uj). Thus Qq{uj) = ~e'^VG{<^)uG{'^), 
where vg{uj) = Y.a^i'^Gai^^) = 4i/gq(w). The last 
equation is a consequence of the symmetry of the MZI. 
In the zero frequency limit vg{uj) reduces to vg{0) = 
2LgI {hvp)- The total charge in region G is Qg{^) — 
e^VG{uj){Va,u - 4ug(w)). 

We generalize now to the case were a voltage is applied 
not only to one of the contacts but also to gate G. The 
gate voltage is labeled Vcit)- In this situation the charge 
in arm G is Qg{,<^) = Cg{ug(uj) — Vg,ui)- Combining 
with our previous result for the charge leads to 



Qg{^) = Cg{ug{^)~Vg^u.) 

= e^yG{i-^){Va.u - 4MG(t^)). 



(20) 



Solving for the internal potential and invoking the defi- 
nition of the injectivity vci}^) (see Eq. |l^) allows us to 
express the internal potential euG{^) through the applied 
voltages: 



ug(^) = 



-iojCGVc^u. + eV(2/t) (1 - e'""'^) Vg, 
-iuoCG + (2e2//i) (1 - e*"^G) 



(21) 



The current in gate G is given by Ig,u = i^QG{^)i where 
—Qci^) is the charge accumulated in the gate. Since, 
with the help of Eq. ( ^l| ) we can express Qg{^) as a func- 
tion of external voltages only, we can calculate the con- 
ductance matrix elements Ggg{'^) = dlG,uj/dVG,uj and 
GGa(^) = dlG,uj I dVa^i^. Notc that the matrix elements 
Gab{^) and Gba{^) vanish since the charge in region G 
is independent of the voltage applied to the gate further 



away from it. (This is a consequence of our assumption 
of forward scattering only at the junctions and of the ab- 
sence of capacitive coupling between the two arms) . For 
later use we here state the result for Ggg{^), which is 



Ggg(^) 



dir 



—iioCG 



dVG,u 1 - 2iuGGRql (1 - e*"-G) • 



(22) 



In Eq. (22) we have introduced the charge-relaxation re- 
sistance Rq = h/{Ae^X-pi the interferometer. The charge 
relaxation resistancellS is a measure of the dissipation 
generated by the relaxation of excess charge on the con- 
ductor into the reservoirs. For a structure with perfect 
channels connected to a reservoir each reservoir channel 
connection contributes with a conductance Gq = 2e^ /h: 
the conductances of different channels add in parallel 
since each channel reservoir connection provides an addi- 
tional path for charge relaxation. For example a ballistic 
wire connected to two reservoirs has a charge relaxation 
resistance Rq = {Gq + Gq)~'^ = h/ie^. For the MZI 
considered here an excess charge in the upper or lower 
branch has the possibility to relax into the four reservoirs 
of the MZI. But at each junction the two connections 
are only open with probability T and R (see Appendix 
B). Thus the two connections act like one perfect chan- 
nel. As a consequence the charge relaxation resistance 
for our MZI is just that of a perfect wire and also given 
by Rq = h/iie^). 
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FIG. 3. The spectrum Eq. (|25| ) is shown as a function of 
the parameter u)t. The full curve corresponds to an interac- 
tion strength g — 0.6, for the middle (dashed) curve g — 0.7 
and for the lowest (dotted) curve g — 0.8. 



In the low-frequency limit we get from Eq. (|2 
Ggg(w) 



-iGfj^.Gi^ 



(23) 



G 7j2^3 
~^q'-'ti,G 



i0' + 



This is in agreement with the result of Blanter et al.E£l 
for a single wire coupled to a gate. We have introduced 
the dimensionless (Luttinger) parameter qg as a measure 
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of the strength of couphng between arm G and gate G. 
If arm and gate are decoupled the interaction parame- 
ter takes the value go = ^ while it goes to zero as the 
strength of coupling is increased. The parameter go is 
related to the capacitance Cg and to the density of states 
Dg = >^g{<-^ = 0) = 2Lg/ (hvp) of the wire throughll3 



9h = 



I + c^Dg/Gg' 



(24) 



The electrochemical capacitancell^ 

c;:}, = Ca' + {e^DGr\ 



Ga,G of arm G is 



The remaining elements of the conductance matrix, 
namely those involving the currents in the contacts of 
the MZI, can be derived from an ac-scattering approach. 
These calculations are presented in Appendix |A[ To 
discuss the influence of potential fluctuations on dc- 
transport we need the spectrum SuoUai^) of these fluc- 
tuations. Since the spectrum of the current fluctuations 
Si^ii^{uj) in region G is related to the real (dissipative) 
part of the element Ggg{<^) of the emittance matrix 
through the fluctuations dissipation theorem Sj^i^{uj) = 
2fcTReG'GG(^) (in the limit hu; < fcT), we get SucUci'^) 
from the relation Su^uG{^)=^iGiGi^)/ (^^^g)'- 



{l-glf{l-cos{LOTG)) 



2(1 - 9gY (1 - cos(wtg)) + 252,(1 - gy)ujTG sin(wTG) + g%{ujTGf 



(25) 



The spectrum Eq. ( [25| ) is shown in Fig. (||) as a function 
of the dimensionless parameter lvtg for different values of 
the interaction parameter gG- Zeroes of SucUci^) occur 
when LOT is a multiple of 27r. This is a consequence of 
our approximation which considers only uniform poten- 
tial fluctuations. 
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FIG. 4. The exact expression for the spectrum Eq. ( |2q ) 
(full line) is compared to the expression Eq. (|2^) (dashed line) 
in the weak coupling limit. Here (? = 0.9 is chosen for the in- 
teraction parameter. 



If the spatial dependence of the potential is taken into 
account, Blanter et al.E3 find that the traversal time is 
renormalized through the interaction (tg — > 5gtg) and 
consequently the zeroes of Su ^(iS) are shifted accord- 
ingly. Instead of the dynamics of single carriers it is plas- 
mons which govern the high frequency dynamics. This 
comparison indicates thus the limitation of our approach: 
Since we start from a one-particle picture our approach 
is most reliable in the case of weak coupling ^g 1- In 
the weak coupling limit [gG — + 1) we expand the spec- 
to the leading order in G^i,^g/Gg — (1 — .9g) 



trum Eq. 
which leads to 



(26) 



The spectrum vanishes in the non-interacting limit gG = 
1 iGf,,G/CG = 0). In Fig. I the full spectrum Eq. (H) 
is compared to the approximate form Eq. (|2^). The 
function sin(wrG/2)/(wrG/2) reflects the ballistic flight 
of carriers through an interval of length L. In the limit 
of strong coupling the potential noise is white and 



SugUg{^) = "ikTRg = 2kTh/{Ae^). 



(27) 



Remarkably in the strong coupling limit for the ballis- 
tic ring considered here the spectrum is universal. The 
only property of the system which enters is the number 
of leads which permit charge relaxation. 

Having found the fluctuation spectra of the internal 
potentials we are now in the position to investigate the 
influence of fluctuations on dc-transport. 



V. ATTENUATION OF THE 
DC-CONDUCTANCE THROUGH COUPLING TO 
A GATE 

We now come to the discussion of the dephasing of the 
coherent part of the dc-conductance in the linear trans- 
port regime due to (equilibrium) charge fluctuations. In 
the last section we have shown that interactions lead to 
effective charge and potential distributions in the arms 
of the ring which in turn give rise to displacement cur- 
rents in the gates and contribute to the ac component 
of the particle currents in the contacts of the ring. The 
zero frequency contribution to the currents in the con- 
tacts remained unchanged. Here we go one step fur- 
ther and will discuss how charge fluctuations act back 
on the dc-conductance. In contrast to the last section 
we will thus only discuss the zero frequency component 
of the conductance. Electron-electron interactions affect 
the coherent part of the dc-conductance only. This can 
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be understood from the well known result0 that interac- 
tions do not change the conductance of a one-dimensional 
wire attached to reservoirs. In the interferometer, when 
interactions are considered, AB oscillations of the con- 
ductance are suppressed. The dc-conductancc matrix 
for the case without interactions is given in appendix 
(see Eq. ([A9|)). Throughout this section we choose 
= fJ'O + e.V and = Ms = = A*o- We will first 
assume that only arm A is coupled to a gate {Cb — > oo). 
The generalization to the case where both arms are cou- 
pled to gates is straight forward and will be discussed at 
the end of this section. 

We will from now on treat the potential as a function 
with certain statistical properties. The potential fluctua- 
tions are characterized by the spectrum Su^u^i^) which 
is defined through 



2nS{uj + uj')SuaUa{^) = (uA(t^)uyi(w'))- 



(28) 



The spectrum was evaluated in Sec. ^. In ad dition the 
potential has a vanishing mean value (C/4(t)) = 0. Here 
(..) denotes the statistical average and the ua{oj) are the 
Fourier components of UA{t) (cf. Eq. (|3|)). The bar is to 
emphasize the distinction between quantum averages (..) 
and statistical averages. 

The quantity of interest to us is the statistically aver- 
aged dc conductance, defined through 



(29) 



A convenient starting point for the calculation of the con- 
ductance the following expression for the current in a 
contact aca, 

t(i) = iy dEdE'e'^^-^"^'/'' (30) 

X [d)^{E)d^{E')-bi{E%{E')\ , 

in terms of the operators aHE) {aa{E)) creating (de- 
stroying) an electron in a state with energy E entering the 
system through contact a and the operators ba(E) and 
ha{E) respectively creating and annihilating an electron 
outgoing at energy E. The operators da{E) and ha{E') 
are related through the scattering matrix (see Fig. (||) 
and Eq. (0)) 



E 



dE'S^p{E,E')di3{E') 



(31) 



As described in Section |l| Eqs. {^,^) and (0) we de- 
termine the scattering matrix elements from the WKB 
solution of the Schrodinger equation for the arm of the 
ring. Doing this we go to the second order in the per- 
turbing potential. This is necessary since due to the as- 
sumption of a vanishing mean value of the statistically 
averaged internal potential there exist no first order cor- 
rections to the ave rag ed conductance Eq. (^9|). Combin- 
ing Eqs. (|3^) and ( ^l| ) with the scattering matrix (|^) and 



statistically averaging leads to the following expression 
for the average conductance: 



(G„0) = ^ / dE(-^^ ME))- 



(32) 



In above equation we have introduced the statistically 
averaged transmission probability 



ME)) 



l±cos(e(S)-27r$/$o) 



(33) 



X {l-(rX,iLA,t))/2h^) 



T sin (eiE) - 27r$/$o) {rAMLA, t))/h 



Here Q = k^.ELA — kB^ELs is a geometric phase, $ 
is the magnetic flux enclosed by the ring and $o is 
the flux quantum. The upper sign is for the pairs 
of indices (a,/3) = (1,3), (2, 4), the lower sign is for 
(a,/3) = (1,4), (2, 3). Furthermore, {T^piE)) and 
{Tpa{E)) are related via the Onsager relations. Ex- 
pressions for r^_i(Lyi,t) and rA,2{LA,t) are given in 
Eqs. (^), (^ and (0). If charge fluctuations are not 
ta ken into a ccount the transmission probability simply 
is {TcAE)) = (1 ± cos {Q{E) - 2Tr$/$o)) /2 (compare 
Eq. (A9)). Interactions thus decrease the amplitude 
of the AB-oscillations and lead to an additional out-of- 
phase contribution. Eq. ( p3| ) can be rewritten in an ap- 
proximate but more convenient form as 



ME)) 



1 r 
2 



(34) 



cos (q{E) - 2^$/$o + {rA,2{LA.t))/h 



Note that also in the presence of interactions cur- 
rent is conserved and the system is gauge invariant. 
This is reflected in the fact that {Gap) = and 
UpW^ = (with JG^ = -eyn ) . Eq. dH ) 
has a rather intuitive interpretation since (TpaiE)) ^ 
(|*^(x,t) + ^'|(a;,t)|2), where 'i>i{x,t) and ^'|(x,t) are 
the WKB wave-functions for the upper and lower arm re- 
spectively at energy E (see Eq. (||)). Note, that Eq. ( |3^ ) 
is strictly correct only to second order in the fluctuating 
potential. 

The transmission probability {Taf3{E)) depends on en- 
ergy E only through the geometric phase Q — kA^sLA — 
kB.sLB ~ —kA.E^L + AELa/ i^vp)- Here it is as- 
sumed that both the difference in length of the two arms 
AL = La — Lb and the difference of the sub-band ener- 
gies in the two branches of the ring AE = Ea — Eb are 
small. The AB oscillations in the conductance are washed 
out completely by thermal averaging when Q(E + kT) — 
Q{E) « 27r. It follows that oscillations may be observed 
at temperatures kT <C Ep{Xp/ AL). Assuming that AL 
is of the order of the Fermi- wavelength Ai? we conclude 
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that kT <C Ep is sufEcient for neglecting the influence 
of thermal averaging on the phase-coherence of the AB- 
oscillations in the conductance. This rather surprising 
result is a consequence of the absence of closed orbits in 
the interferometer. It is generally expected that meso- 
scopic phase coherence is destroyed at temperatures of 
the order of the Thouless energy Et — Ef{Xf/L) ^ Ep- 
If closed orbits were considered the transmission proba- 
bility would be a function not only of kpAL but also 
of kpL. However, in our simple model we can write 
(Gap) w {e^/h){TapiEF)) for kT < Ep (but still in 
the limit huj <C kT) when assuming that AL ^ Xp. 




FIG. 5. The ratio of traversal time and dephasing time 
t/t^ versus the coupling parameter g. The upper (full) line 
shows the numerically evaluated exact result (cf. Eq. (^)), 
while the dashed line is the weak coupling approximation 
Eq. (^). In the calculations we used the parameters of 
Ref. [8], namely La ~ 1.5/im, kp ~ 1.5-10*m~^. Furthermore 
we assumed T — IK and for the effective mass m* of the elec- 
tron we used the value for GaAs, m* = 0.067me. With these 
parameters the Fermi velocity is vp ~ hkp/m* ~ 2.6 • 10^ 
m/s. 

We can now proceed to further evaluate Eq. (|3^) . Us- 
ing rA,i{LA,t) as given in Eq. (0) and invoking the defi- 
nition of the spectrum Eq. ( p8| ) we obtain 

1~2 — /r A 2 f sin^(wr/2) 

J ZTT LjJ 

The fluctuation spectrum SjjauA^) has been calculated 
in Sec. The integral in Eq. ( p5| ) can be done analyti- 
cally in the two limiting cases of very strong [g — > 0) and 
weak {g — > 1) coupling between arm and gate. Since the 
spectrum SuaUa{^) (Eq. (p5|)) is a function of u) only 
through LOT, it can be seen from Eq. (^5|) that we can 
then write 

{r\ALA,t))/{21i?) = taIt^, (36) 

where ta is the traversal time and is a function of 
temperature and the coupling parameter only. Eq. ( |3^ ) 
defines the dephasing time r^. It is expressed through 



the phase-shift r^^i(L^, t) acquired by a WKB wavefunc- 
tion in the presence of a time dependent potential (rel- 
ative to the case without potential) and quantifies the 
strength of the suppression of the AB oscillations in the 
dc-conductance (see Eq. (|4|)). In the limit where gate 
and ring are weakly coupled we use the approximate spec- 
trum Eq. (|2^) to evaluate Eq. (^5|). The dephasing rate 
~ '''(^^ found to be 

r*-(J)f<'-»J)' = (i)1f'™'- 

The dephasing rate = t^^ is linear in tempera- 
ture. Very recently experimental results were reported 
by Hansen et aLQ on the temperature dependence of de- 
coherence of AB-oscillations in ballistic rings. In these 
measurements a dephasing rate linear in temperature 
was found. The dephasing rate (Eq. (^) also de- 
pends on the coupling parameter gA- Dephasing goes 
to zero when gate and ring are completely decoupled 
{gA — > 1). We can similarly determine the dephasing 
rate in the strong coupling limit. We know that in this 
case the potential noise is white and the spectrum is 
given by SuaUa{'^) = '^kTRq. The dephasing rate is 
= Rq{e'^ /h^)kT. Due to the more complicated form 
of the fluctuation spectrum SuaUa ('^) i'^ the intermediate 
parameter range for g we can't give a simple analytical 
expression for the dephasing rate. However, the integral 
in Eq. ( ^5| ) can be performed numerically. In Fig. (|^) 
t/t^ is plotted versus the interaction parameter gA over 
the full parameter range. 

Returning to the weak coupling limit we further eval- 
uate rA,2{LAit) and find 

5Q = {rAMLA,t))/h = - (^^) . (38) 

To be consistent with our previous approximations this 
term should be neglected since it is of the order kT/Ep. 
Still it is interesting to compare the size of the two cor- 
rections due to scattering from the internal potential. We 
find that (r/r^) /(56 ~ kpL 1 which imphes that tak- 
ing scattering into account to the first order in the po- 
tential is a surprisingly good approximation. Combining 
the information gathered so far allows us to rewrite the 
transmission matrix elements Eq. (|^) in the more con- 
venient form 

(TapiE)) - ^ (l ± e-"/"* cos {e{E) - 27r$/$o))) , 

(39) 

where the energy dependent part of Q{E) is of the or- 
der kT/Ep in the limit where the two arms have similar 
length AL ^ Xp. The theory developed so far can be 
readily adapted to the case where both arms are coupled 
to gates {gc ^ Q,G = A, B) by making the replacement 
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TB 



T(I,,A T^^B 



(40) 



in Eq. (pQI). The simple result Eq. (|40|) is an immedi- 
ate consequence of the fact that potential fluctuations in 
the two arms are uncorrelated ( SuaUb — SubUa — 0)- 
This can either be seen from the corresponding matrix 
elements of the admittance matrix (cf. Appendix ^ or 
by directly calculating potential correlations, as is done 
in the low frequency limit in Appendix H. 



VI. CONCLUSIONS 

In this work we have examined dephasing due to 
electron-electron interactions in a simple Mach-Zehnder 
interferometer (MZI). Without interactions the MZI ex- 
hibits only forward scattering. (However, screening will 
generate displacement currents at all contacts in response 
to a carrier entering the conductor). We have shown how 
a scattering matrix approach can be used to calculate the 
effect of charge fluctuations on the conductance. We have 
first determined the internal potentials and their statis- 
tical properties in a non-perturbative way. Subsequently 
we calculated corrections to the dc-conductance up to the 
second order in the effective internal potentials. In the 
expression for the averaged dc-conductance a dephasing 
time T0 occurs in a natural way. It is a measure of the 
strength of the attenuation of the AB-oscillations as a 
function of temperature and coupling strength between 
ring and gates. The dephasing rate = ^ is found to 
be linear in temperature and to depend on the Luttinger 
coupling parameter g through cc (1 — g^)^, at least 
in the weak coupling limit. Alternatively, it depends on 
the ratio of the electrochemical capacitance and the 
geometrical gate capacitance C like (C^/C)^. In terms 
of the Coulomb energy Ec = e^/2C needed to charge 
the wire and the density of states (inverse level spacing) 
D = 2L/hvF this ratio is (1 + l/{2DEc)y^. Such a de- 
pendence on Ec cannot be obtained from a Golden rule 
argument in which the coupling between the ring and the 
gate is treated perturbatively. Such a treatment would 
lead to a dephasing rate proportional to E^- A dephasing 
rate proportional to E'^ is obtained only in the (unrealis- 
tic) limit that the level spacing far exceeds the Coulomb 
energy. On the other hand we found that the evaluation 
of the phase accumulated during traversal of the conduc- 
tor is surprisingly well described just be the first order 
perturbation theory in the fluctuating potential. 

Recently the temperature dependence of AB- 
oscillations in ballistic fjngs was investigated experi- 
mentally by Casse et al.El. Since both thermal averag- 
ing and dephasing of the electronic wavefunctions lead 
to a decrease in the visibility of the AB-oscillations as 
temperature is increased a separation of these different 
effects is of interest. Such an analysis of experimental 
data was carried out by Hansen et aha. They find that 
the dephasing rate is linear in temperature in agreement 



with our work. The dephasing length (Z^ = vft) we 
have calculated can be of the or|der of the dephasing 
length observed in this experimentu when the coupling g 
is taken to be strong enough. A more detailed compari- 
son of our result to the experiment is difficult, since the 
setup of Ref. || is different from MZI presented here. In 
the experimental setup a top-gate is used to cover the 
(two-terminal) AB-ring and no side-gates are used. 

We note here only as an aside that a linear temper- 
ature dependence was also observed in experiments on 
chaotic cavitieaH3. The theory presented here, i.e. the 
spatially uniform fluctuations of the potential in the in- 
terior of the cavity, will Also give rise to a linear in tem- 
perature dephasing rateO. It is also interesting to note 
that our dephasing time shows features similar to the 
inelastic scattering time for a ballistic one-dimensional 
wire. The inelastic scattering time of Ref. 24 is inversely 
proportional to temperature and can be written as a sim- 
ple function of the Luttinger liquid parameter measuring 
the strength of electron-electron interactions. Whether 
the inelastic time of Ref. ^ is in fact also the dephasing 
time remains to be clarifled. 

Our discussion has emphasized the close connection 
between the ac-conductance of a mesoscopic sample and 
dephasing. We have given the entire ac-conductance ma- 
trix for the model system investigated here. A current 
and charge conserving ac-conductance theory requires a 
self-consistent approach to determine the internal poten- 
tials and requires the evaluation of the displacement cur- 
rents. 

The displacement currents at the gates can in principle 
be measured. Nevertheless the thermal potential fluctu- 
ations inJjhe arms of the ring do not act as a which path 
detcctorEj. In fact the dephasing rate increases with de- 
creasing capacitance and is maximal if C = 0. In this 
limit there are no displacement currents at the gates. 
The absence of which path detection is reflected in the 
fact that the charge correlations in the two arms of the 
ring vanish in the equilibrium state of the ring. 

The work presented here can be extended in many 
directions. Multi-channels systems and systems with 
backscattering can be considered. The role of the ex- 
ternal circuit can be examined. We hope therefore that 
the work reported here stimulates further experimental 
and theoretical investigations. 
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APPENDIX A: ADMITTANCE MATRIX 



ico) = {L (^)) = /° {co) + II (^) + (^). (A4) 



We here give the admittance matrix Gik{uj) ~ 
dh.ui / dVk^ui {I, k — 1,2, 3, 4, A, B) for a symmetric {La = 
Lb = L and thus t^i = = r) Mach-Zehnder interfer- 
ometer (MZI). The resuhs for the asymmetric case are 
similar but notationaUy more cumbersome. As in the 
rest of the paper we consider the hmit Tilij <C kT and 
huj,kT ^ Ep- We have shown in Sec. |^ how the ad- 
mittance matrix elements that relate the displacement 
currents in the gates to voltages applied at a gate or a 
contact can be calculated. It remains to derive the ma- 
trix elements that relate currents in the contacts to ex- 
ternal voltages. To this end we employ the ac-scattering 
approach following Ref. O. We here give a slightly for- 



19 and 



malized derivation of the results found in Ref 
generalize the results to accommodate a system like the 
MZI containing several regions described by different in- 
ternal potentialscJ (the two arms in the case of the MZI). 
For recent related work we refer the reader to Refs. ^ 28 
A time-domain version of the ac-scattering approach was 
recently introduced in Ref. to the investigation of 
charge pumping in open quantum dots. 

We consider the situation where a time-dependent volt- 
age Vfj, (t) is applied to contact {^) of the system. We start 
from the current operator in the form Eq. (|30|). Fourier 
transforming gives 



(Al) 



- bi{E)K{E + huj) . 



To take into account scattering from internal potential 
fluctuations we use Eq. (|3l|) and write 



= e dE 



..(EK(S + M-/m (A2) 



J2 SlpiE, Ei)S^^{E + hu:, L;2)4 (^1)0^(^2) 



Our next step is to average this expression quantum me- 
chanically. Doing this it has to be taken into account, 
that the distribution of charge carriers coming in from 
reservoir is modified due to the time-dependent volt- 
age apphed to this contact (see Eq.^. Since we consider 
only the linear response we expand the scattering matrix 
to first order in the internal potentials ug{u)). We write 

Sap{E', E) = 27,n5{E - E')S^^l{E) + sl'^{E', E), (A3) 

where S'^^^ (E) is the scattering matrix of the ideal ballis- 
tic system (see below Eq. (|])) and S^^p{E' , E) is linearly 
proportional to the perturbing potential. Substitu ting 
Eqs. ( [Tsl ) and (A3) into the current operator Eq. (A2) 
and taking the quantum average we see that average cur- 
rent in contact a may be written in the form 



The first term /^(cj) in Eq. (A4) is the dc contribution 
/°(c^)=e<5(M (A5) 



xJ2 I dE 



0-\s'fl{E)\ 



In the case of interest to us here, fp{E) — f{E) and 
/^(w) = due to the unitarity of the scattering matrix 

S^^^{E). The current 

.f{E)- f{E + hu;)^ 



flLJ 



-V0{u), 



can be understood as the part of the total current di- 
rectly injected into contact a due to the oscillations of 
the external potential V^{t) (see Ref. |9|). The third con- 
tribution to the total current is the response to 
the internal potential distribution (compare also Ref. |l^) 

X {f{E)- f{E + hLj)). 

We now want to apply the theory developed so far to 
the calculation of the dynamic conductance of the MZI. 
For this example the full scattering matrix is given in 
Eq. ( I31I) . Inelastic transitions are absorbed in the scat- 
tering matrices of the arms Sa{E',E) and Sb{E',E). 
From Eqs. (@), and (|ll|) we know that to first order 
in the potential 



Sg{E + TiLU, E) = 27r(5(u;)e*'=^-^<^ 



(A8) 



gifcEiG!^£M (i_e*' 



where we used kE+hu ~ fc^; + cj/u. Expressions for 
the matrix elements S^^\E',E) for the inte rfer ometer 
are now easily derived from Eqs. ( pT| ) and (AS). For 
the calculation of the admittance matrix it is further- 
more important t o no te, th at in the limit of interest here 
{huj ^ kT) Eqs. (A6) and (A7) can be considerably sim- 
plified. The Fermi functions appearing in these equations 
are expanded to linear order in huj/kT. Since in addition 
the products of sca ttering matrix elements contained in 
Eqs. (A6) and (A7) do not depend on the energy E but 
only on the energy difference E' — E = fno for the scat- 
tering matrix used here, the energy integrations can be 
performed. 

Combining the scattering matrix elements defined in 
Eq. (0) with the expressions for the currents in the 
gates Ig{^) = itoQcito) and the currents in the con- 
tacts Ia{uj) = I^{uj) + I^{oj) (see Eqs. (M) and (0)) we 
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can now calculate all elements of the conductance matrix. 
We here consider the case of a perfectly symmetric (La = 
Lb) interferometer and give the result to first order in oj. 
Expanding the admittance matrix in the low frequency 
limit wc can write Gikioj) — Gff. — iujEik + o{uj'^). Ex- 
plicitly, the zeroth order term is 



/Go 







^G^ 


o\ 





Go 


^G^ 


^G^ 







^G^ 


Go 








^G^ 


^G^ 





Go 




















Vo 











0/ 



Go 



where we have introduced 



(A9) 



Go = e^h, 

G| -eV(2/i)(l±cos(27r$/$o)). 



(AlO) 
(All) 



Note that in the dc-limit there are no currents in the 
gates. The first order term Eik is called the emittance 
matrix. It is of the form 



E : 



/ 




-E 


Et 


E^ 


-Ea 


—Eb 
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-E 


E^ 


Et 


-Ea 


—Eb 






Et 
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-E 


-E 


-Ea 


—Eb 










-E 


-E 


-Ea 


—Eb 
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-Ea 


-Ea 
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-Eb 


—Eb 


-Eb 


—Eb 





^Eb 


I 



(A12) 



The entries of the omittance matrix are defined 
through 



E ■ 

Eg = 
E^ = 



g (G^,,a/Ga + C^^b/Gb) t , 
±E^ + E + EaI2 + Eb/2, 



Eq, = — rcos($). 
The electrochemical capacitance G~^ 



(A13) 

(A14) 
(A15) 

(A16) 



G 



is defined as 
where Dq — 2Lq / {hvp) is 



the density of states per unit length. The charge relax- 
ation resistance Rq is given by Rq — /i/(4e^) and the 
traversal time is tq = Lq/vp. Current conservation im- 
plies Gikioj) = while as a consequence of gauge 
invariance Gik{i^) = 0. Similar sum rules hold for 
every coefficient in the expansion of Gik{uj) as a function 
of u; (e.g. Elk = 0, ^Li Eik - 0). 



APPENDIX B: CHARGE-CHARGE 
CORRELATIONS 

The spectra of charge fluctuations in the gates or, 
equivalently, the arms of the interferometer, as well as 
correlations between the charges in the two gates, can be 



calculated directly from the knowledge of the scattering 
matrix Sap without first calculating the dynamic conduc- 
tance as we have done in Sec. IV. This approach which 



is particularly convenient in the low frequency limit was 
introduced in Ref. ^ for a mesoscopic structure coupled 
to a single gate and generalized to the case of coupling 
to more than one gate in Ref. pl| . In this section we 
apply the approach of Refs. 30 sTto calculate the zero 



frequency limit of the fluctuation spectra of the charges 
in the gates of the MZI. In contrast to the rest of this pa- 
per we calculate both equilibrium and non-equilibrium 
fluctuations. From Ref. ^ it is known that in equilib- 
rium, to leading order in frequency, the charge-charge 
correlations are given by 



^QgC 



2kTC^.GC^.HR. 



GH 



(Bl) 



where G and H specify the gates in the system. Here 
G 



Cq -\- {e Dq)~^ is the electrochemical capaci- 
tance of gate G, Dg being the density of states. Fur- 
thermore i?^^ is the generalized charge relaxation resis- 
tance to be introduced below. For G = H Eq. (Bl) gives 
the spectrum of charge fluctuations in gate G while for 
G H Eq. (Bl) gives the equilibrium charge correlations 
between gates G and H. With a small voltage applied 
to one contact of the system the fluctuation spectra to 
leading order in the applied voltage are Ref. |l| 



qV 



2\eV\G^,GC^,HR^". 



(B2) 



The generalized charge relaxation resistance Rf^ and 
the corresponding quantity in the driven case, the Schot- 
tky resistance Ry^ can be expressed through the (off- 
diagonal) elements of a generalized Wigner Smith time- 
delay matrix: 



GH 



r: 



pGH _ 



2? DgDh ' 
, T nG jjH* 



2e2 



DgD 



H 



(B3) 
(B4) 



The density of states Dg in region G is the sum of the di- 
agonal matrix elements of the Wigner Smith time-delay 
matrix, Dg — J2a Eaa- The elements of the time-delay 
matrix can conveniently be found from the scattering ma- 
trix via the relationEa 



Da^{E) 



1 



2TTi ^ 

13 



^SUE)'^ 



(E) 



dE 



(B5) 



The scattering matrix for the MZI is given in Eq. (|^). We 
here only need the matrix sj^p [E) which is found from 
Eq. (§) by replacing Sg{E' , E) in Eq. (g) by Sg{E' , E) = 



d{E' 



E)Sg 



(E) with S^q\e) = expiikpLG) ( here 
we simply have Qg = ^Lg)- To calculate the time- 
delay matrix in region G we replace the energy deriva- 
tive by a derivative with regard to the local potential 
d/dE —d/dUG- Since the scattering matrix depends 
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on energy onlyj-through the phase factors 8g 
is easy to showE£l that 



= -ttDc 



dS 



01 



de 



G 



kLc it 



(B6) 



We can now use Eqs. (B3) to (B6) to calculate the gen- 
erahzed charge relaxation and Schottky resistances: 

' (B7) 



Kg 



R 



BB 



R 



AB 



R 



BA 



4e2' 
0, 



nAA 



R 



AB 



— jyBB 

— ^1/ 



= = - 



16e2' 
h 



(B8) 
(B9) 

(BIO) 



Combining these equations with Eqs. (Bl) and ( p32| ) we 
get charge-charge correlations for the equihbrium and 
out-of-cquihbrium situations. It is interesting to note 
that in equihbrium the charges in the two gates are un- 
correlated. This is a consequence of the absence of closed 
electronic orbits in the ring and the fact that we have not 
introduced a Coulomb coupling between the two branches 
of the ring. For the same reason correlations are inde- 
pendent of the magnetic field. This implies, that despite 
the fact that AB-oscillations are observed in the currents 
measured at one of the contacts, the interior of the ring 
behaves like a classical system. In contrast, the charge 
fluctuations generated by the shot noise are correlated. 
Like the equilibrium charge fluctuations they are for the 
geometry investigated here independent of the AB-flux. 
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